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Analytic structure of the gluon and quark propagators in Landau gauge 
QCD * 

R. Alkofer, Institute for Theoretical Physics, University of Tubingen, D-72076 Tiibingen, Germany, 
W. Detmold, Department of Physics, University of Washington, Box 351560, Seattle WA 98195, USA, 
C. S. Fischer, IPPP, University of Durham, Durham DHl 3LE, U.K., 

P. Maris, U. of Pittsburgh, Dept. of Physics and Astronomy, 100 Allen HaU, Pittsburgh PA 15260, USA 

In Landau gauge QCD the infrared behavior of the propagator of transverse gluons can be analytically deter- 
mined to be a power law from Dyson-Schwinger equations. This propagator clearly shows positivity violation, 
indicating the absence of the transverse gluons from the physical spectrum, i.e. gluon confinement. A simple 
analytic structure for the gluon propagator is proposed capturing all important features. We provide arguments 
that the Landau gauge quark propagator possesses a singularity on the real timelike axis. For this propagator we 
find a positive definite Schwinger function. 



The standard model of particle physics consists 
of gauge field theories. These have been postu- 
lated on the basis of symmetries and their ele- 
mentary excitations do not reflect the observed 
particle spectrum. In the quantum formulation 
of these theories, especially in Poincare-covariant 
gauges, an intricate problem is posed by the sepa- 
ration of physical and unphysical degrees of free- 
dom. 

In Quantum Electrodynamics (QED) in linear 
covariant gauges, the electromagnetic field can 
be decomposed into transverse, longitudinal and 
timelike photons, however, only transverse polar- 
izations are observed. From a purely mathemati- 
cal point of view, this can be understood from the 
representations of the Poincare group for massless 
states: massless particles have only two possible 
polarizations 1 . The apparent contradiction is 
resolved by the fact that timelike and longitudi- 
nal photons cancel exactly in the S'- matrix 1^. 
In this context we emphasize that the states of 
quantum gauge field theories in covariant gauges 
necessarily constitute an indefinite metric space. 
In covariant gauge QED one thus has to sacrifice 
the principle of positivity of the representation 
space. 



* Summary of a talk given at several occasions; to be pub- 
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In Quantum Chromodynamics (QCD) in lin- 
ear covariant gauges, the cancellation of unphys- 
ical degrees of freedom in the iS-matrix is sub- 
stantially complicated by the self-interaction of 
the gauge fields and by the ghost fields that are 
necessarily present in the quantum formulation 
of these theories . To order Og in perturbation 
theory, one obtains amplitudes for the scattering 
of two transverse gluons into one transverse and 
one longitudinal gluon. However, at the same or- 
der, a ghost loop appears and cancels the various 
gluon loops, and scattering of transverse to lon- 
gitudinal gluons does not occur. It is possible to 
prove this cancellation to all orders in perturba- 
tion theory on the basis of the BRS 4 symme- 
try of the covariantly gauge fixed theory. This 
symmetry can be represented by gauge transfor- 
mations with the ghost field as a parameter. The 
ghost fields, being scalar, anti-commuting fields, 
are necessarily in the unphysical part of the rep- 
resentation space. In covariant gauge QCD, the 
physical (and thus positive definite) part of the 
state space is conjectured to be the set of BRS 
singlets 

Gauge fixed QGD is invariant under transfor- 
mations related to the ghost number, and employ- 
ing such transformations, one can show that BRS 
non-singlets occur in quartets Such a BRS- 
quartet consists of two parent and two daugh- 
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ter states of respectively opposite ghost numbers. 
The latter states are BRS-exact and thus BRS- 
closed because the BRS transformation is nilpo- 
tent. The BRS daughters are orthogonal to all 
other states in the positive definite subspace and 
therefore do not contribute to physical S'-niatrix 
elements. The parent states belong to the indef- 
inite metric part of the representation space and 
are expected to violate positivity. 

The Kugo-Ojima confinement scenario '6J5j de- 
scribes a mechanism by which the positive semi- 
definite, physical state space contains only color- 
less states. Colored states are not BRS singlets 
and therefore do not appear in S'-matrix elements: 
they are confined. Thus an investigation of (non- 
) positivity of transverse gluons and quarks allows 
us to understand confinement via the BRS quar- 
tet mechanism in more detail. 

Here we present analytic properties of the gluon 
and quark propagators in Landau gauge QCD 
as they result from non-perturbative calculations 
(more details can be found in Ref. [?])• We con- 
firm previous results I8l9j on positivity violation 
for the gluon propagator. We also provide a pa- 
rameterization of the gluon propagator that is an- 
alytic throughout the complex plane except 
on the real timelike axis and which decreases to 
zero in every direction of the complex plane. 
Such behavior satisfies the usual axioms of lo- 
cal quantum field theory JOI (except positivity). 
For the quark propagator, we analyze several gen- 
eral constraints, lattice data JI], and solutions of 
the coupled quark-gluon-ghost Dyson-Schwinger 
equations (DSEs) None of these contradict 

positivity of the quark propagator. 

Within the framework of a Euclidean quantum 
field theory (used in the following), positivity is 
formulated in terms of the Osterwalder-Schrader 
axiom of reflection positivity jlHI . In the special 
case of a two point correlation function, A(a; — 
y), the condition of reflection positivity can be 
written as 



d^x d^y fix, ~xo) A{x - y) f{y, yo) > , (1) 



where f(x,xo) is a complex valued test function 
with support in {{x, xq) : a;o > 0}. After a three- 
dimensional Fourier transformation, this condi- 



tion can be given in terms of the Schwinger func- 
tion 

A(i) = ^ J dpcos{tp)a{p^) >0, (2) 

where crip^) is a scalar function extracted from 
the corresponding two-point function. 

The elementary two-point functions of QCD 
are the ghost, gluon, and quark propagators. 
In Landau gauge these renormalized momentum- 
space propagators Dg{p), D^^{p), and S{p) can 
be generically written as 

G(p2) 



Dg{p) = — 



pi 



S{p) = 



PiJ,Pu\ Z{p'^) 
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(3) 
(4) 



-i-^A{p^) + B{p^) 

= i^a^{p^)+as{p^). (5) 

Note that the quark propagator S{p) is decom- 
posed into scalar and vector parts, Osip^) and 
a^ip^)- Violation of reflection positivity can be 
studied by calculating Eq. Q with (Jg{p^) = 
Z{jP')/p'^ for the transverse gluons and (Js{p^) and 
<Jy{p'^) for the quarks. (The ghost propagator vi- 
olates reflection positivity by the way ghosts are 
introduced in Faddeev-Popov quantization.) 

The coupled set of DSEs (for recent reviews 
see e.g. the ghost and gluon propagators 

can be solved analytically for p^ 0+ JHl ■ One 
finds simple power laws for the gluon and ghost 
dressing functions 

Z{p') ^ {p^ f-, G{p') ^ {p^r\ (6) 

such that the product Z{k^)G^{P) ^ as{k^) 
goes to a constant in the infrared: there is an in- 
frared fixed point for the running coupling. The 
value of the exponent k is in the range 0.5 < 
K < 0.7, depending on details of the truncation 
of the set of DSEs ;12.l6ll7ll8j . Here w e use 
a self-consistent truncation scheme |12ll(j| which 
neglects the effects of the four-gluon interaction 
and employs ansdtze for the ghost-gluon, quark- 
gluon, and three-gluon vertices such that two im- 
portant constraints are fulfilled: the running cou- 
pling, asip'^), is independent of the renormal- 
ization point, and the anomalous dimensions of 
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Figure 1. DSE [El and lattice results for the 
gluon renormalization function Z{p'^) and the fit 
(jHll are shown. 



the propagators are reproduced at one-loop level 
for large momenta. In this particular truncation 
K = (93 - yi20T)/98 « 0.595 [TTIT^ is an ir- 
rational number, and as(0) « 2.972. Infrared 
dominance of the gauge fixing part of the QCD 
action '19' implies infrared dominance of ghosts in 
the DSEs which, in turn, can be used to show |17| 
that as (0) depends only weakly on the dressing of 
the ghost-gluon vertex and not at all on other ver- 
tex functions. Recently, the same infrared behav- 
ior of the propagators has been found using the 
method of Exact Renormalization Group Equa- 
tions EOl- 

For quenched QCD, the gluon propagator, as 
it results from numerical solution of the coupled 
DSEs for the gluon and ghost propagators, agrees 
very well with recent lattice data see Fig. 
The unquenched DSE gluon propagator is signifi- 
cantly suppressed in the intermediate momentum 
region, where screening effects of qq pairs becomes 
important. For both Nf = and Nj = 3 the in- 
frared behavior is given by Eq. ©. 

The corresponding running coupling can be ac- 
curately represented by ^21 



with Pq = [llNc-2Nf)/2>. Note that the Landau 
pole at spacelike = Aqcq is canceled, c./. Ref. 
|22| . The expression ((TJ is analytic in the complex 

plane except on the real timelike axis {p^ < 0) 
where the logarithm produces a cut. 

Since the infrared exponent, k, is an irrational 
number, one knows that the corresponding gluon 
propagator possesses a cut on the negative real 
axis as well. It is possible to fit the DSE solution 
for the gluon propagator without introducing fur- 
ther singularities. The fit to the gluon renormal- 
ization function with the form 

with w — 2.5 and Aqcd = 510 MeV, is in good 
agreement with the DSE solution. Here w is 
a normalization parameter, and 7 = (— 13A'c -t- 
4Nf)/{22Nc - 4Nf) is the one-loop value for the 
anomalous dimension of the gluon propagator. 
The discontinuity across the cut in Zfit(p^) van- 
ishes for ^ 0~, diverges to -l-oo at = ~^qcd 
on both sides and drops to zero for p^ — > —00. 
Additional parameterizations have been explored 
in Ref. 7,. 

The corresponding Schwinger function, Ag{t), 
based on the fit, Eq. JHJ, is compared to the DSE 
solution in Fig. 13 To enable a logarithmic scale, 
the absolute value is displayed. Ag{t) has a zero 
for i « 5 GeV~^ « 1 fm and is negative for larger 
Euclidean times. I.e. we clearly observe positiv- 
ity violations in the DSE gluon propagator and 
the agreement of the numerical Schwinger func- 
tion with the Fourier transformed fit is also ex- 
cellent. The crucial property for positivity viola- 
tion in the gluon propagator is that it vanishes 
for ^ 0+. This can be seen from the relation 
= a-g{p^ = 0) = J d'^x D{x) which implies that 
the gluon propagator in coordinate space, D{x), 
is trivially zero or necessarily contains positive as 
well as negative contributions. 
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Figure 2. Shown are the results for the absolute 
value of the gluon Schwingcr function from the 
DSE solution, the fit 0, and the infrared part of 
this fit. 



In expression the overall magnitude, w, is 
arbitrary because of renormalization properties 
(it is determined via the choice of the renormal- 
ization scale). The infrared exponent, k, is deter- 
mined from the infrared analysis of the DSEs, and 
the one-loop value is used for the gluon anomalous 
dimension, 7. Thus, we have found a parameter- 
ization of the gluon propagator which has effec- 
tively only one parameter, the scale Aqco. This 
and the relatively simple analytic structure gives 
us confidence that we have succeeded in uncover- 
ing the important features of the Landau gauge 
gluon propagator. 

We now turn to the analytic structure of the 
quark propagator. In the DSE studies of Ref. 
|12| it has been assumed that the non-Abelian 
part of the quark-gluon vertex can be factored out 
from the tensor structure. This structure is then 
given by the Curtis-Pennington vertex \2'^\ which, 
by construction, is multiplicatively renormaliz- 
able and satisfies the Abelian Ward-Takahashi 
identity of QED. In this particular case one has 
a term 
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Figure 3. Results for the quark Schwingcr func- 
tions using different approximations to the quark- 
gluon vertex. 



in the vertex. For the present discussion it is im- 
portant to note that the exact quark-gluon vertex 
(which satisfies the more complicated Slavnov- 
Taylor identity of QCD), and future improved 
approximations to it, will almost certainly also 
exhibit such a quark-gluon coupling proportional 
to the sum of quark momenta. This coupling, be- 
ing effectively scalar, is not invariant under chiral 
transformations in contrast to the leading 7^ part 
of the vertex. However, the term in Eq. (^J only 
appears once chiral symmetry is broken dynami- 
cally and is thus consistent with the chiral Ward 
identities. Its existence provides a significant self- 
consistent enhancement of dynamical chiral sym- 
metry breaking in the quark DSE, which is neces- 
sary to produce an acceptable value of the chiral 
condensate ^21 ■ Quite independently of the form 
of the gluon propagator, the resulting quark prop- 
agator respects positivity if the term in Eq. 
is included in the quark-gluon vertex [7], as evi- 
denced by the Schwingcr functions in Fig. O (solid 
and dashed curves). From this figure we also see 
that positivity is violated if only the bare quark- 
gluon vertex is used (dotted curve), in agreement 
with previous work employing this approxima- 
tion 1211. 
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A single real pole on the negative momentum 
axis results in a pure exponential decay of the cor- 
responding Schwinger function. This is not what 
is observed in (t) for small times where there is 
some curvature (see Fig.|31)- To give an acceptable 
reproduction of (t) , analytic parameterizations 
of the quark propagator must contain additional 
sub-dominant structure. By Fourier transforming 
combinations of functions with cuts and poles at 
timelike momenta we have determined three sim- 
ple sources for such curvature. The dominant sin- 
gularity may be accompanied by additional real 
singularities at larger mass scales, or by complex 
conjugate singularities with a larger real part of 
the mass, or it may be the starting point of a 
branch cut on the negative real momentum axis 
Of course such a list is by no means exhaus- 
tive. 

In Ref. we have used various general re- 
quirements, lattice data ^J, and the DSE so- 
lutions |12| to constrain these various paramet- 
ric forms. Such analysis leads to an important 
conclusion: in all parameterizations, the domi- 
nant singularity must occur on (or very close to) 
the real, timelike half-axis at a scale m 350 to 
500 MeV, the lower estimate coming when lattice 
data are used for the fit, the upper estimate for 
fits to the DSE solutions. This scale may relate 
to a constituent quark mass. 

Whilst this result is robust, our current meth- 
ods are not able to accurately determine the na- 
ture of the dominant singularity or reliably con- 
strain the additional sub-dominant structures. In 
particular, the regular infrared behavior of quark 
propagator found in the DSE and lattice calcu- 
lations makes determining its detailed analytic 
properties very difficult without explicitly prob- 
ing the timelike momentum half-plane. A more 
conclusive analysis of the structure of the quark 
propagator would result from direct solution of 
the DSEs over an appropriate region of the com- 
plex momentum plane. 

To summarize: We have proposed the rela- 
tively simple function, Eq. © , to describe the full 
(non-perturbative) Landau gauge gluon renor- 
malization function for all complex values of mo- 
mentum. The corresponding gluon propagator 
agrees well with DSE solutions and lattice data 



for spacelike momenta, it is positivity violating, 
and it is analytic everywhere except for a cut on 
the negative real half-axis. Thus it implies 
that in Landau gauge QCD, the confinement of 
transverse gluons is related to the violation of 
Osterwalder-Schrader reflection positivity. We 
have also provided evidence that the Schwinger 
functions related to the quark propagator are pos- 
itive definite, and consequently quark confine- 
ment is not manifest at the level of the propa- 
gator. A number of relatively simple parameter- 
izations have been suggested for this propagator 
in terms of real and complex conjugate poles and 
branch cuts. In all cases, the dominant singular- 
ity is real and occurs at a scale 350 to 500 MeV. 
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